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Abstract 

We have considered a hot strange star matter, just after the cohapse of a supernova, as a 
composition of strange, up and down quarks to calculate the bulk properties of this system at finite 
temperature with the density dependent bag constant. To parameterize the density dependent bag 
constant, we use our results for the lowest order constrained variational (LOCV) calculations of 
asymmetric nuclear matter. Our calculations for the structure properties of the strange star at 
different temperatures indicate that its maximum mass decreases by increasing the temperature. 
We have also compared our results with those of a fixed value of the bag constant. It can be seen 
that the density dependent bag constant leads to higher values of the maximum mass and radius 
for the strange star. 



* Corresponding author 

'I' E-Mail: bordbar@physics.susc.ac.ir 

^ Permanent address 



1 



Keywords: Strange star, equation of state, structure, density dependent bag constant 



INTRODUCTION 



Strange stars are those which are built mainly from self bound quark matter. The surface 
density of strange star is equal to the density of strange quark matter at zero pressure 
(~ 10^^ g/cm?), which is fourteen orders of magnitude greater than the surface density of a 
normal neutron star. The central density of these stars is about five times greater than the 



surface density 



l|-|3| . The existence of strange stars which are made of strange quark matter 



was first proposed by Itoh [4| even before the full developments of QCD. Later Bodmer jsl 
discussed the fate of an astronomical object collapsing to such a state of matter. In 1970s, 
after the formulation of QCD, the perturbative calculations of the equation of state of the 
strange quark matter was developed, but the region of validity of these calculations was 
restricted to very high densities j6|. The existence of strange stars was also discussed by 
Witten 0]. He conjectured that a first order QCD phase transition in the early universe could 
concentrate most of the quark excess in dense quark nuggets. He suggested that the true 
state of matter was strange quark matter. Based on theoretical works of Witten on cosmic 
separation of phases, the transition temperature is approximately 100 MeV , an acceptable 
QCD temperature [7|. Witten proposal was that the strange quark matter composed of light 
quarks is more stable than nuclei, therefore strange quark matter can be considered as the 
ground state of matter. The strange quark matter would be the bulk quark matter phase 
consisting of almost equal numbers of up, down and strange quarks plus a small number of 
electrons to ensure the charge neutrality. A typical electron fraction is less than 10~^ and it 
decreases from the surface to the center of strange star [l-3]. Strange quark matter would 
have a lower charge to baryon ratio compared to the nuclear matter and can show itself in 



the form of strange stars 



7H10| 



Just after the collapse of a supernova, a hot strange star may be formed. A strange star 
may be also formed from a neutron star and is denser than the neutron star. If sufficient 
additional matter is added to a strange star, it will collapse into a black hole. Neutron 
stars with masses of 1.5 — I.8M0 with rapid spins are theoretically the best candidates for 
conversion to the strange stars. An extrapolation based on this indicates that up to two 
quark-novae occur in the observable universe each day. Besides, recent Chandra observations 
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indicate that objects RX J185635-3754 and 3C58 may be bare strange stars 

In this article, we consider a hot strange star born just after the collapse of a supernova. 
Here we ignore the effects of the presence of electrons, and consider a strange star purely 
made up of the quark matter consisting of the up, down and strange quarks. The energy 
of quark matter is calculated at finite temperature, and then its equation of state is de- 
rived. Finally using the equation of state of quark matter, the structure of strange star at 
different temperatures is computed by integrating the Tolman-Oppenheimer-Volkoff (TOV) 
equations. 

II. CALCULATION OF QUARK MATTER EQUATION OF STATE 
A. Density Dependent Bag Constant 

Different models have been used for deriving the equation of state of quark matter. 
Therefore there is a great variety of the equations of state for this system. The model which 
we use is the MIT bag model which was developed to take into account the non perturbative 
effects of quark confinement by introducing the bag constant. In this model, the energy per 
volume for the quark matter is equal to the kinetic energy of the free quarks plus a bag 
constant (B) 12|. The bag constant B can be interpreted as the difference between the 
energy densities of the noninteracting quarks and the interacting ones. Dynamically it acts 
as a pressure that keeps the quark gas in constant density and potential. This constant 
is shown to have different values which are 55 and 90 in the initial MIT bag model. 
Since the density of strange quark matter increases from surface to the core of the strange 
star, it is more appropriate to use a density dependent bag constant rather than a fixed bag 
constant. 

According to the analysis of the experimental data obtained at CERN, the quark-hadron 



transition takes p 
(156 MeVJm-^) 



for B 



15 



ace at about seven times the normal nuclear matter energy density 



13 



, 1^. Recently, a density dependent form has been also considered 
18|. The density dependence of B is highly model dependent. In this article, the 
density dependence of B will be parameterized, and we make the asymptotic value of B 
approach a finite value Boo |l8| , 
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Bin) = i3oo + (Bo - i3oo)e~^("/"«)'. (1) 

The parameter Bq = B{n = 0) has constant value which is assumed to be Bq = 400 in 
this work, and 7 is the numerical parameter which is usually equal to hq ^ 0.17 fm~^, the 
normal nuclear matter density. Boo depends only on the free parameter Bq. We know that 
the value of the bag constant {B) should be compatible with experimental data. The exper- 
imental results at CERN-SPS confirms a proton fraction Xpt = 0.4 (data is from experiment 
on accelerated Pb nuclei) Q- Therefore, in order to evaluate ^Boo, we use the equation 
of state of the asymmetric nuclear matter. The calculations regarding this can be found in 
the next section. 



B. Computation of Boo using the asymmetric nuclear matter calculations 

We use the equation of state of the asymmetric nuclear matter to calculate Boo- For 
calculating the equation of state of asymmetric nuclear matter, we employ the lowest order 



constrained variational 



energy as follows 19l427|. 



XOCV) many-body method based on the cluster expansion of the 



The asymmetric nuclear matter is defined as a system consisting of Z protons (pt) and 
neutrons (nt) with the total number density n = Upt + Unt and proton fraction Xpt = 

where Upt and n„t are the number densities of protons and neutrons, respectively. For this 

system, we consider a trial wave function as follows, 

^ = Fcj), (2) 

where is the slater determinant of the single-particle wave functions and F is the A-body 
correlation operator {A = Z + N) which is taken to be 

F=5n/(u) (3) 

i>j 

and iS is a symmetrizing operator. For the asymmetric nuclear matter, the energy per 
nucleon up to the two-body term in the cluster expansion is 



1 < 


H 


ijj > 


A Kip 


ip > 
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The one-body energy, Ei, is 

1=1 ki 

where labels 1 and 2 are used for proton and neutron respectively, and ki is the momentum 
of particle i. The two-body energy, i?2, is 

^2 = ^E<*J|V(12)|2J-J2>, (6) 

where 

V(12) = -|^[/(12), [V?2, /(12)]] + /(12)1^(12)/(12) • (7) 

In the above equation, /(12) and V{12) are the two-body correlation and nucleon-nucleon 
potential, respectively. In our calculations, we use UVu + TNI nucleon-nucleon potential 
2^ . Now, we minimize the two-body energy with respect to the variations in the correlation 
functions subject to the normalization constraint. From the minimization of the two-body 
energy, we obtain a set of differential equations. We can calculate the correlation functions 
by numerically solving these differential equations. Using these correlation functions, the 
two-body energy is obtained and then we can compute the energy of asymmetric nuclear 



matter. The procedure of these calculations has been fully discussed in reference 201]. 



As it was mentioned in the previous section, the experimental results at CERN-SPS 



confirms a proton fraction Xpt = 0.4 13|, [iSj, therefore to compute Boo-, we proceed in the 
following manner: 

• Firstly, we use our results of the previous section for the asymmetric nuclear matter 
characterized by a proton fraction Xpt = 0.4. By assuming that the hadron-quark 
transition takes place at the energy density equal to 1100Mey/m~^ Q, Q, we find 
that the baryonic density of the nuclear matter is = 0.98/m~^ (transition density). 
At densities lower than this value the energy density of the quark matter is higher than 
that of the nuclear matter. With increasing the baryonic density these two energy 
densities become equal at the transition density, and above this value the nuclear 
matter energy density remains always higher. 

• Secondly, we determine Boo = 8-99 by putting the energy density of the quark 
matter and that of the nuclear matter equal to each other. 
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C. Calculations for the energy of quark matter at finite temperature 



To calculate the energy of quark matter, we need to know the density of quarks in terms 
of the baryonic density. We do this by considering two conditions of beta equilibrium and 
charge neutrality. This leads to the following relations 

(J-d = f^u — /^e, (8) 

f^s f^u f^ei (9) 
fJ'S = fJ'd, (10) 

2/3nu- l/Sris - l/3nd- He = 0, (11) 
where /ij and rii are the chemical potential and the number density of particle i, respectively. 



As mentioned, we consider the system as pure quark matter (rig = ) [8|, l29l-l3lj . Thus 
according to relation flTT]) . we have 

n^ = l/2{n, + nd). (12) 

The chemical potential, /ij, at any adopted values of the temperature (T) and the number 
density (n^) is determined by applying the following constraint, 

g roc 

n. 



271^ Jo 

where 



roo 

/ f{ni,k,T)edk, (13) 

•J 



/(n„ T) ^^p^^^y ^ n^k^c^y/^ - fi,)} + 1 ^^^^ 



is the Fermi-Dirac distribution function 32|. In the above equation, (3 = and g is the 
degeneracy number of the system. 

As it is previously mentioned, we consider the total energy of the quark matter as the 
sum of the kinetic energy of the free quarks and the bag constant (B). Therefore, the total 
energy per volume of the quark matter (Stot) can be obtained using the following relation, 

Stot = £u + Sd + £s + B, (15) 

where Si is the kinetic energy per volume of particle i, 

^^ = 7^ r + h^k^c^f'fin,. k, T)edk. (16) 



After calculating the energy, we can determine the other thermodynamic properties of 
the system. The entropy of the quark matter {Stot) can be derived as follows 



Si 



'tot — 



(17) 



where Si is the entropy of particle i, 



Siirii.T) 



--kB [f{n„k,T)\n{f{n,,k,T)) 

TT^ JO 

+ (1 - /(n„ k, T)) ln(l - /(n„ k, T))]k''dk. 



(18) 



The Helmholtz free energy per volume (J-") is given by 



■tot ~ 



TSi 



'tot- 



(19) 



The entropy per particle of the quark matter as a function of the baryonic density for 
two cases of the constant and density dependent B at different temperatures are plotted 
in Figs. [1] and [21 For a fixed temperature, we see that the entropy per particle decreases 
by increasing the baryonic density and for all relevant densities, it is seen that the entropy 
increases by increasing the temperature. 

In Figs. [3] and IH the free energy per volume of the quark matter versus the baryonic 
density for two cases of the constant and density dependent B are presented at different 
temperatures. We can see that the free energy of the quark matter has positive values for 
all densities and temperatures. For all densities, it is seen that the free energy decreases by 
increasing the temperature. 

To obtain the structure of the strange star, the equation of state of the quark matter is 
needed. For deriving the equation of state, the following equation is used. 



where P is the pressure. The pressure of the quark matter versus the baryonic density for 
two cases of the constant and density dependent B are plotted in Figs. |5] and |6l It is seen 
that by increasing both density and temperature, the pressure increases. These figures show 
that for each temperature, the pressure becomes zero at a specific value of the density. We 
see that the density corresponding to zero pressure increases by decreasing the temperature. 



(20) 
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III. STRUCTURE OF STRANGE STAR 



Compact objects like white dwarfs, neutron stars and strange stars have hmiting masses 
(maximum mass) and with a mass more than the hmitting value, the hydrostatic stability 
of the star is impossible. For obtaining the maximum mass of the strange star, we use the 



Tolman-Oppenheimer-Volkoff (TOV) equations |29| . 



dP_ G[g(r) + ^][m(r) + ^^] 

^ = ATTr'^Sir). (22) 
dr 

By using the equation of state found in the previous section, we integrate the TOV equations 
to calculate the structure of the strange star {2^. The results of this calculation are given 
in the following figures and tables. 

Figs. [7] and [8] show the gravitational mass versus the central energy density at different 
values of temperature for two cases of the constant and density dependent B. For each 
value of the temperature, these figures show that the gravitational mass increases rapidly 
by increasing the energy density and finally reaches to a limiting value at higher energy 
densities. It is seen that the limiting value of the gravitational mass increases by decreasing 
temperature. Comparing Figs. [7] and El one concludes that at all temperatures, for the 
density dependent bag constant, the rate of increasing mass with increasing the central 
density, at lower values of the central densities, is substantially higher than that of the case 
for fixed bag constant, especially at zero temperature. In Figs. [9]and[l0l we have plotted 
the radius of strange star versus the central energy density for both B = 90^^^ and density 
dependent B at different temperatures. From Figs. ITHTOl it can be seen that at each central 
density, both mass and the corresponding radius increase by decreasing the temperature. 
The gravitational mass of strange star is also plotted as a function of the radius for the 
constant and density dependent B in Figs. [11] and [121 It is seen that for all temperatures, 
the gravitational mass of strange star increases by increasing the radius and it approaches a 
limiting value (maximum mass). Figs. [TT] and [T^ show that by decreasing the temperature, 
the limiting values of mass and the corresponding radius both increase. 

In Tables [T] and [Til the maximum mass and the corresponding radius and central energy 
density of the strange star at different temperatures for two cases of the constant and density 
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dependent B are given. It is shown that by decreasing the temperature, the maximum 
mass of strange star increases. This behavior is also seen for the radius of strange star 
versus the temperature. Meanwhile, the central energy density decreases by decreasing the 
temperature. By comparing Tables [T] and HIl we can see that for all temperatures, the 
maximum mass and the corresponding radius calculated with the constant B are less than 
those calculated with the density dependent B. 

IV. SUMMARY AND CONCLUSION 

We have considered a pure quark matter for the strange star to calculate the structure 
properties of this object at finite temperature. For this purpose, some thermodynamic 
properties of the quark matter such as the entropy, free energy and the equation of state 
have been computed using the constant and density dependent bag constant {B). It was 
shown that the free energy of the quark matter decreases by increasing the temperature while 
the entropy of this system increases by increasing the temperature. It was indicated that 
by increasing the temperature, the equation of state of the quark matter becomes stiffer. 
We have calculated the gravitational mass of the strange star by numerically integrating 
the Tolman-Oppenheimer-Volkoff (TOV) equations. Our results show that the gravitational 
mass of the strange star increases by increasing the central energy density. It was shown 
that this gravitational mass reaches a limiting value (maximum mass) at higher values of the 
central energy density. We have found that the maximum mass of the strange star decreases 
by increasing the temperature. It was also shown that the maximum mass and radius of the 
strange star in the case of density dependent B are higher than those in the case of constant 
B. 
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TABLE I: Maximum mass {Mmax) in solar mass unit (Mq), and the corresponding radius (R) and 
central energy density (Sc) of the strange star at different temperatures (T) for ;B = 90 



T (MeV) MmaxiMQ) R (km) £,{10^'^^) 






1.354 


7.698 


38.24 


30 


1.228 


7.073 


47.54 


70 


1.101 


6.416 


60.60 


80 


1.039 


6.142 


63.65 



TABLE II: As Table U but for the density dependent B. 



T (MeV) MmaxiM^) 


R (km) 







1.676 


8.761 


39.11 


30 


1.341 


7.442 


48.47 


70 


1.181 


6.768 


61.56 


80 


1.122 


6.567 


64.21 
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FIG. 1: The entropy per particle of the quark matter versus the baryonic density at different 
temperatures for B = 90 
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FIG. 2: As Figured] but for the density dependent B. 
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FIG. 3: The free energy per volume of the quark matter versus the baryonic density at different 
temperatures for 5 = 90 ^j^- 
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FIG. 4: As Figure [3] but for the density dependent B. 
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FIG. 5: The pressure of the quark matter as a function of the baryonic density at different tem- 



peratures for ^ = 90 
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FIG. 6: As Figure [5] but for the density dependent B. 
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FIG. 7: The gravitational mass of the strange star as a function of the central energy density at 



different temperatures for ;B = 90 
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FIG. 8: As Figure [7] but for the density dependent B. 
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FIG. 9: The radius of the strange star as a function of the central energy density at different 
temperatures for B = 90 ^j^- 
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FIG. 10: As Figure E] but for the density dependent B. 
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FIG. 11: The gravitational mass of the strange star as a function of the radius at different tem- 
peratures for ^ = 90 
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FIG. 12: As Figure [TT] but for the density dependent B. 



21 



